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Abstract 

We present a non-perturbative study of the phase diagram of SU(2) Yang-Mills theory 
in a five-dimensional spacetime with a compact extra dimension. The non-renormalizable 
theory is regularized on an anisotropic lattice and investigated through numerical simulations 
in a regime characterized by a hierarchy between the scale of low-energy physics, the inverse 
compactification radius, and the cutoff scale. We map out the structure of the phase diagram 
and the pattern of lines corresponding to fixed values of the ratio between the mass of the 
fifth component of the gauge field and the non-perturbative mass gap of the four-dimensional 
modes. We discuss different limits of the model, and comment on the implications of our 
findings. 



1 Introduction 



Gauge theories denned in a spacetime with compact extra dimensions offer a very interesting plat- 
form for model building and phenomenology beyond the Standard Model, and provide appealing 
theoretical scenarios for Grand Unification and electroweak symmetry breaking, for the fermion 
hierarchy problem, and for the strong CP problem [TMl4j. 

Although a gauge theory in d > 4 is generally non-renormalizablc, and thus cannot be considered 
as a fundamental theory, it can be interpreted as a low-energy effective description of an underlying 
theory The effective theory is consequently a description valid only for energies lower than a cutoff 
A, which is the scale where the microscopic details of the underlying theory become relevant and 
the effective description breaks down. In particular, for a compact extra dimension of size L5, 
this effective theory description is expected to be sensible only if there is at least a modest scale 
separation between the inverse cutoff A -1 and the compactification length, L5A ^> 1. 

Assuming that the relevant correlation lengths are much longer than the inverse cutoff A -1 , 
the actual ultra-violet completion of the theory can be traded for any simpler regularization, 
introducing only finite systematic errors. In particular, choosing the lattice regularization, one can 
obtain robust non-perturbative predictions for the low-energy sector of the effective theory from 
numerical simulations. In this work, we investigate a simple model of this type: non-Abelian SU(2) 
Yang-Mills (YM) theory in a flat (4+l)d Euclidean spacetime. 

The properties of 5d YM models on the lattice have been studied in several works, including 
Rcfs. [15H36] . In fact, our approach builds on the evidence from previous numerical studies, showing 
that SU(2) and SU(3) YM theories defined on an isotropic 5d space exhibit a Coulomb phase in the 
weak (bare) coupling regime [T5HT8] . and that upon compactification of one of the directions, the 
theory undergoes dimensional reduction for any compactification radius [19p20| . It is the existence 
of the Coulomb phase and dimensional reduction that guarantees that the relevant correlation 
lengths are large compared to the cutoff and thus justifies the use of the lattice regularization. 

In our model, the dimensional reduction allows one to describe the low-energy sector of the 
(4+l)d theory using 4d YM theory, which has a non-perturbatively generated string tension pro- 
portional to 



where (75 denotes the gauge coupling of the (4+l)d theory. Measuring all quantities in units of 
the string tension (instead of the cutoff) allows one to build a consistent definition of the lattice 
spacing and of lines of constant low-energy physics, in terms of 4d quantities. 

This paper is structured as follows: after defining the lattice formulation of the model in Sec- 
tion [21 we study the phase diagram of the theory in Sec. [3J where we also discuss some qualitative 
analogies of the model with compact U(l) lattice gauge theory in (3+l)d. Our investigation of the 
properties of the dimensionally reduced phase is reported in Sec. [4J while in Sec.[5]we comment on 
possible strategies of removing the cutoff and why they must eventually fail. Finally, we summarize 
our findings and present an outlook on future research directions in Sec. [6] Some technical details 
of our calculation are discussed in the Appendix. 



As a prototype for an extra-dimensional gauge model, we study SU(2) Yang-Mills theory in (4+l)d 
Euclidean spacetime, where one of the dimensions is taken to be periodic with a finite extent L5 
and the gauge fields arc taken to obey periodic boundary conditions along the compact direction. 
The model is formally defined by the path integral 



where the M, N indices run from 1 to 5. Since the bare coupling g\ has the dimensions of a 
length, this model is non-renormalizable, and the continuum action in Eq. ^ defines it only up 
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2 The model 
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to a rcgularization, i.e. as a theory with a cutoff A. The model is thus parametrized by two 
dimcnsionless ratios, namely L$A and A. 

To study the model at non-zero coupling values, we regularize it on an anisotropic lattice with 
a lattice spacing a in the four usual directions, letting the lattice spacing as in the fifth direction 
take an independent value. Furthermore, let N5 denote the extent of the compact direction in 
units of 05, such that L5 = (I5N5. We define the model on the lattice through the Wilson action 



~ x,l<M<N<A 



x,KM<4 



l--ReTrP MS (a:) 



(3) 



with 



A = (4) 
In the weak coupling limit, the anisotropy factor 7 reduces to the ratio of the lattice spacings 

lim — =7. (5) 

The non-renormalizability of the model implies that the limit in which both a and 05 are taken 
to zero leads to divergences which cannot be absorbed in coefficients of an action with a finite 
number of terms, and thus the cutoff cannot be removed. However, it is possible to take the lattice 
spacing in the compact direction to zero, 05 — s- 0, while keeping ^5/7 = N5 w L$/a and the 
lattice spacing a in the other four directions fixed0 In this work, we present results obtained from 
extrapolations to this limit, which offers two advantages: 

• it reduces the number of parameters characterizing the model from three (N$, 7, and ,85) 
down to only two (N$ and /?s) like in the continuum theory; 

• it suppresses cutoff effects in systems where the L^/a ratio is not very large. Practically, this 
is very important in our numerical calculations, since the 4d correlation length of the system 
strongly depends on this ratio. 

After taking the continuum limit in the fifth dimension, the 4d lattice spacing a cannot be taken 
to zero, but one can still approach the continuum limit of the theory by tuning the correlation 
length of the lattice model £ to be much larger than the lattice spacing a: in such a limit, the 
low-energy sector of the model becomes insensitive to the lattice discretization artifacts, and the 
continuum 4d symmetries are restored. 

We have investigated the phase diagram and various correlation lengths of the model non- 
perturbatively via Monte Carlo simulations, calculating two-point correlation functions of Polyakov 
loops along the extra dimension, as well as along the four ordinary directions (torelons). In par- 
ticular, the torelon mass E t , the string tension er, and the energy of the first Kaluza-Klein excita- 
tion Ekk can be extracted from the two-point correlation function of zero transverse-momentum 
torelons P x (r) along one of the space-like directions of size L: for asymptotically large separations 
r one has 

(Pi(r)P a (0)) * Aie~ EtV + A 2 e~ E ^ r , (6) 

with 



L 6 ' V 1 \L 



-5 



Note that, besides the term proportional to the string tension, E t also includes a Liischer-term 
correction [371135], which is the leading contribution coming from the massless fluctuations of the 
torelon along the dt = 2 transverse dimensions (see also Ref. [33] for a discussion). 



1 This leads to a quantum-mechanical system defined on a 4d lattice with finite spacing a, and is not in conflict 
with non-renormalizability issues. 
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Figure 1: The effect of the discretization due to finite N 5 on the pseudo-critical coupling f3 c (N 5 /^). 
The phase transition breaks the center symmetry of the Polyakov loop along the extra dimension, 
such that (P5) vanishes in the confining phase but has a non-zero expectation value in the di- 
mcnsionally reduced phase. The pseudo-critical coupling is defined through the maximum of the 
susceptibility X 5 = (|A| 2 > - (|A|) 2 - 



3 Phase diagram 

The phase structure of 5d SU(2) Yang-Mills theory was first studied non-perturbatively in Ref. [T5] 
for an isotropic action on a lattice with same extents in all directions. Unless a dimensionful scale L5 
is introduced by compactification, the theory has only one parameter, namely the lattice coupling 
05. In this case, the phase diagram consists of a confining phase in the strong-coupling regime at 
small /?5 values, and a Coulomb- like phase at weak coupling, i.e. at large ,85. The confining phase 
is separated from the Coulomb phase by a transition, which (based on numerical evidence from 
simulations on a 4 5 lattice) Ref. [15] located at /3| = 1.65. This is identified as a first-order bulk 
phase transition, revealed by a discontinuity in the plaquette expectation value and not related 
to any breaking of a global symmetry. As the confining phase is not analytically connected to 
the weak-coupling regime, we consider it merely as a lattice artifact, having little to do with the 
continuum theory. 

Since in the absence of a compactification the weak-coupling phase is characterized by an infinite 
correlation length, we expect that the compactification of the fifth dimension leads to spontaneous 
breaking of the center symmetry and to dimensional reduction, for any finite compactification 
radius. For this reason we denote this phase, characterized by a non-zero Polyakov loop expectation 
value ((P5) 7^ 0), as the "dimcnsionally reduced phase". Our numerical results indeed confirm the 
existence of two phases also for the anisotropic model defined by Eq. ([3]), for all compactification 
lengths N$. 

In order to map out the phase diagram, we measured the pseudo-critical coupling for fixed N*, 
and for various N$ extrapolating to the N§ — >• 00 limit, depicted in Fig. d]. The resulting phase 
diagram in this limit is displayed in Fig. [2] where the two phases are separated by a second order 
transition line, instead of a first order transition as in the isotropic case. The order and universality 
class of the phase transition were established by a finite-size scaling study discussed in more detail 
in the Appendix. At the transition, the correlation length of Polyakov loops in the fifth direction 
diverges, and, as expected from general arguments |40j . we found the transition to be in the trivial 
4d Ising universality class. 

Finally, it is interesting to point out that some features of the (4+ l)d YM system are analogous 
to the more familiar case of an Abelian, U(l) gauge theory defined on a (3 + l)d lattice, whose 
properties are well understood [4TT443j . 

It is well known that compact U(l) lattice gauge theory in 4d admits a dual formulation, see, 
for instance, Ref. [44] and references therein, in terms of topologically non-trivial excitations which 
can be localized in elementary lattice cubes |45j and interpreted as magnetic monopoles. Numerical 

2 The discretization errors are of order I/N5 as they arise from dimension six operators, some of which have 
coefficients proportional to aa^. 
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Figure 2: The phase diagram of 5d SU(2) lattice gauge theory with a compact extra-dimension, in 
the N$ 00 limit: it consists of a confining and of a dimcnsionally reduced phase, separated by a 
second-order phase transition for which the Polyakov loop along the compact extra dimension serves 
as the order parameter. The solid red curve, drawn to guide the eye, goes through the origin with 
divergingjlerivative, as expected from the strong-coupling expansion, and tends to an asymptotic 
value as N5 — > 00. The confining phase is analytically disconnected from the weak-coupling, large 
/?5 regime, and is an artifact of the lattice regularization. 

investigations [35] have revealed that such monopoles are responsible for the phase structure of 
the theory. In the weak-coupling regime ((3 = 1/e 2 3> 1) at zero temperature, the monopoles have 
a large mass m(/3), and their worldlines (which form closed loops, due to current conservation) 
are typically short, since the probability of a loop of size L is proportional to exp[— m(j3)L] <C 1. 
Such loops do not affect the long-distance physics: in the presence of external electric charges, the 
system exhibits Coulomb behavior (like its counterpart in the continuum limit, namely a photon 
gas), and the effect of the monopoles is simply to renormalize the electric charge. 

As /3 is reduced and the bare coupling is made stronger, however, the monopoles become lighter, 
and their mass vanishes at a weakly first-order transition at j3 c = 1.0111331(21) [47] — see also 
Rcf. [48] . At stronger coupling, monopoles condense and the system is in a confining phase, 
consistently with the strong coupling expansion of the lattice theory. Thus, the flux of the magnetic 
field coming from a monopole can be seen as disordering space-like Wilson loops. Therefore, at weak 
coupling, monopole loops of finite size can only induce a perimeter-like decay for the expectation 
values of asymptotically large Wilson loops. On the contrary, at strong coupling the monopole 
loops percolate, and disorder the Wilson loops on all scales, driving an exponential area-law decay 
and confining behavior. 

Consider now the compact U(l) lattice gauge theory in a system with a compact fourth direction 
of size L4: in the Coulomb phase at weak coupling, the probability of monopole loops winding 
around the compact direction is proportional to exp[—m{^)L^\. Such loops can disorder spatial 
Wilson loops on all length scales, and thus induce an area-law decay, because the magnetic flux is no 
longer neutralized by that of a neighboring anti-monopolc. As a consequence, spatial Wilson loops 
will exhibit area law behavior, with a string tension proportional to the monopole loop density: 
a 2 a cx exp[— m(/3)L4]. Such behavior has indeed been observed numerically in Ref. [49] : the 
Coulomb phase gives way to a 3d confining phase for any finite L4, where the global U(l) symmetry 
of the Polyakov loop is spontaneously broken. Moreover, by increasing L4 while keeping (3 fixed, one 
can increase the 3d correlation length exponentially: this achieves dimensional reduction and yields 
the continuum limit of the 3d Abelian theory, which is confining for all values of the coupling [50)151] . 

All this is completely analogous to the behavior of YM theory in (4 + 1) dimensions in the 
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Figure 3: Two-point correlation function of zero-momentum torelons as a function of their sep- 
aration, from lattices of sizes 6 2 x L x 24 x 2 with L = 4, . . . , 8, /3 = 1.6575, and 7 = 1. 
The four-dimensional Liischer correction has been subtracted from the correlators such that 
C(r) = \og[(P x (r)P x (0)}]/L— irrj (3L 2 ). The collapse of the curves at large distances indicates that 
the long distance properties of the model are described by a four-dimensional confining theory. At 
short distances the data also reveal the contribution from the first non-static Kaluza-Klein mode. 



weak-coupling phase, although here we have not been able to investigate the objects playing a 
role analogous to that of the magnetic monopolcs as their gauge invariant formulation is yet to be 
found. 



4 Dimensional reduction 

In this Section, we discuss the properties of the dimensionally reduced phase of our model, com- 
paring the pcrturbative predictions which are expected to hold in the weak (bare) coupling limit 
/3s — > 00 and our numerical results from lattice simulations. 

Perturbatively, at large /3s the theory has three mass scales, in analogy with the finite-temperature 
QCD setup in 4d [52VI55] . Firstly, observe that the non-static Kaluza-Klein modes can be described 
by a tower of 4d modes with masses 

2tt , . 

m KK = —n, n = l,2,.... (8) 
^5 

Secondly, at tree-level, the static Kaluza-Klein modes are all massless, but one-loop corrections 
give the gauge field in the fifth direction a mass proportional to 

m 5 cc - w . (9) 

Finally, while the static modes of the gauge field components in the four ordinary directions remain 
massless to all perturbativc orders they acquire a non-pcrturbative mass due to confinement in four 
dimensions. 

In the low-energy sector, namely at distances Ax 3> £5, the system undergoes dimensional 
reduction, and the effective low-energy action describing the static Kaluza-Klein modes takes the 
form of a continuum 4d YM theory coupled to an adjoint scalar field A5 

S eS = J d 4 x Qlx F%, + Tr [D^A 5 ] 2 + m 2 5 Ti A% + ATr A^j , (10) 
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Figure 4: A sample of our results for the string tension era 2 , plotted as a function of /?5. The figure 
refers to data obtained from simulations at a fixed value of N$ = A5/7 = 2, for which Eq. (fTS)) 
predicts a linear relation between \n(aa 2 ) and ,$5. The solid blue line on the right-hand side of the 
plot corresponds the results extrapolated to the N$ — s- 00 limit, and its slope is consistent with the 
perturbative expectation. 



where the effects of the higher order terms are suppressed by powers of the scale ratio between the 
compactification scale (L^/2n) and the typical scale of the effective theory m^ 1 . 

The 4d effective theory given by Eq. (|10[) is rcnormalizable, and is characterized by the param- 
eters g4,([A), m|(/i) and A(/i), where /i denotes the renormalization scale of the effective theory. 
At weak coupling, it is natural to impose a perturbative matching of the parameters and fields of 
the effective theory to the original theory, within the validity domain of the theory described by 
Eq. (fT0|) . Note that Aj$ and g\ get their leading-order contribution at tree- level 



A% = ,/L B A%, 9t=gt/L 5 . (11) 

On the other hand, the leading contributions to m| and A arise only at one-loop level, and 
can be readily evaluated by computing the effective potential using standard techniques. On an 
anisotropic Euclidean lattice in d = 5 dimensions, the one-loop quantum fluctuations induce an 
effective potential for a classical background gauge field parametrized by A$(z) = iTq{z)a^ / '(£555) 
of the form 



y(L) 

►off, 1- 



M . .<=* (^"jf^^M"! 1 ^--^], (12) 

from which the one-loop perturbative result for m| follows: 

ml = \ (^3) K^Lopfe) ■ (13) 

^ 5 ' q=0 

As the 4d components of the gauge field A^ remain massless to all orders in perturbation theory, 
at length scales Ax ^> m^ 1 the A§ field can be integrated out completely, resulting in an effective 
theory which is just a 4d Yang-Mills model 

S' eS = fd*x(±TtFt), (14) 



eff - / " * \ Tj >" 
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Figure 5: The coefficient A describing the exponential decay of the string tension, as a function 
of the fivc-dimcnsional bare coupling era 2 ~ cxp(— Aj3) as a function of N$. The green triangle 
denotes the one- loop pcrturbative prediction as given by Eq. (|T6|) . For lattices coarse in the fifth 
dimension, the coefficient differs by a factor of two, but approaching the continuum limit in the 
fifth dimension improves the agreement significantly. 




Figure 6: Simulation results for the Higgs screening mass E§ extracted from the exponential decay 
of the zero- momentum correlator of the Polyakov loop in the fifth direction, Eq. (|18|) . versus the 
string tension in units of the lattice spacing. The screening mass E§ is expressed in units of the 
square root of the 4d string tension, and is given by = 1m§ to leading order in perturbation 
theory The figure shows our data at fixed N§ — N$/j — 2 and for different N$ values. The 
solid curve is obtained by extrapolation to the N§ — > oo limit. While the data show qualitative 
agreement with the one- loop perturbative prediction (thin dashed line), the quantitative details 
differ due to moderate values of the bare coupling (3. 
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with a gauge coupling <? 4 (/i). 

As before, the coupling of the effective theory can be perturbativcly matched to 54 in the regime 
where both theories are expected to give a good description. Although the correlation length £ 
is infinite to all orders in perturbation theory, the 4d Yang-Mills theory described by Eq. (TT4")) is 
known to be confining for all non-zero values of g' A , since it non-perturbatively develops a mass 
gap, and thus a finite correlation length £ = 1/ ' %fa inversely proportional to the square root of the 
associated 4d string tension a. At leading order, asymptotic freedom requires the dependence of 
the correlation length on the renormalized gauge coupling g'± (/x) to be of the form 



£ = - exp 
H 



1 



(15) 



where bo = ll/(247r 2 ) is the coefficient of the leading term in the /3-function, and c is a constant 
to be determined non-perturbatively. This is the origin of Eq. (fT]). Correspondingly, the large- 
distance decay of Wilson loops normal to the compactified direction is expected to be characterized 
by a string tension 



1 



1 



bog' A \L 



1 



a 2 iV| 



exp 



(16) 



Our simulations show that the scale separation and the dimensional reduction carry over also 
to non-zero couplings. In Fig. [3J we plot zero transverse momentum torclon correlators for var- 
ious torelon lengths L as a function of separation. The four-dimensional Luscher-term has been 
subtracted from the correlators, such that the collapse of the curves at large separations (absent 
if one applies a five-dimensional correction) indicates that the long distance properties of the five- 
dimensional theory are described by a four-dimensional confining theory, as anticipated by the weak 
coupling analysis. While at large distances the correlator is well described the four-dimensional 
theory, at short distances hints of the first non-static Kaluza-Klein mode can be seen. Indeed, 
a two-exponential fit to the data, with the second mass being that appropriate for the first non- 
static mode (given by Eq. ([?])) captures the short distance behavior; the extrapolation of the fitting 
function all the way to zero distance reproduces the data reasonably well. 

Having established the confining nature of the long-distance sector we move on to test the 
effective theory prediction of the behavior of the string tension as a function of the parameters of 
the five-dimensional theory. In Fig. [4] we display a collection of string tensions measured from the 
torclon correlators for several values of N5, keeping N$ fixed and varying the bare coupling /3. For 
all values of N$, the data points form a straight line in a semi- logarithmic plot 



exp(-A0), 



(17) 



as anticipated by Eq. (fl6j). However, for small values of N§ the slope, A, is significantly larger 
than expected from the pcrturbative prediction. Nevertheless, taking the continuum limit in the 
fifth direction brings us significantly closer to the one-loop result as depicted in Fig. [SJ While the 
bare coupling /3 was not that small in the simulations (~ 1.6 — 2) the good agreement between 
perturbation theory and simulations is perhaps not completely unexpected as the coefficient in 
the exponential of Eq. (|16j) is protected from a two- loop correction. The evaluation of the g2(L^) 
to two loops would only induce an overall shift in the amplitude c and thus the first non-trivial 
correction enters at the three-loop level. 

We also measured the correlators of Polyakov loops in the fifth direction from the same data set, 
whose exponential decay 



(H{r)Pm) 



exp(-E 5 r), 



(18) 



is given to leading order in perturbation theory by £5 = 1m§. The results in units of string tension 
are shown Fig. El where it can be seen that the measured data are qualitatively described by the 
leading order perturbation theory, but with sizable quantitative differences. 
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5 Towards the continuum 



Let us now discuss the lines of "constant low-energy physics" in the dimensionally-reduced phase 
of the lattice model, meaning the lines corresponding to a fixed value of the ratio m^j^fa of the 
A§ field mass over the typical energy scale \fa characterizing the 4d physics. In particular, in the 
following we discuss such lines in the phase diagram obtained after taking the continuum limit in 
the fifth direction, N 5 ~ > 00, as a function of fa and N 5 , as sketched in Fig. [7j 

Starting from very large /3g values, leading-order perturbation theory predicts that the lines of 
fixed m^/yja tend asymptotically to hyperboles of equation foN^ = const (in particular, m 5 /y/a 
increases when (3sN 5 takes larger and larger values). As one moves along one such line towards 
smaller values of /3§, the ultraviolet cutoff measured in units of the 4d string tension, era 2 , decreases 
monotonically, and corrections which depend on a/g 2 make the line deviate from its asymptotic 
hyperbolic behavior. Following the line further upwards in the (/3s, -/V5) plane, it eventually ap- 
proaches the phase transition line, and bends to follow it — without crossing it, as long as the latter 
is of second order (as this would imply m^j \fa = 0). 

In Fig. [7] we also display the pattern of lines parametrized by aa 2 = const, i.e. by a fixed value 
of the 4d lattice spacing, as dashed blue lines. The value of a decreases as one moves upwards 
along a vertical straight line, or to the right along a horizontal line. 

Moving along a line of fixed m^/y/a from very large to smaller (3$ values, the lattice spacing a 
decreases. This is not in contradiction with the intuitive expectation that the lattice spacing should 
become smaller when /?5 is increased towards +00, if 7V5 is held fixed. By taking the /?5 — > 00 
limit at fixed N$ (which is described by the green horizontal arrow in Fig. [7]), one crosses lines 
corresponding to larger and larger vn^j \fa values. Note that this corresponds to moving along the 
extrapolated violet curve appearing in Fig. [51 going towards the vertical axis era 2 = 0. In such a 
limit, the A§ static modes (and the Kaluza-Klein tower) become infinitely heavy compared to the 
4d energy scale, and decouple from the low-energy sector of the theory. Thus, in this limit the 
model simply reduces to the ordinary continuum 4d YM theory, which is renormalizable, and for 
which there is no obstruction in removing the ultraviolet cutoff a — > 0. 

More generally, this "ordinary" 4d continuum limit can be reached by considering any trajectory 
that crosses lines corresponding to increasing values of the m^/^/a ratio, while increasing the 4d 
correlation length £ in units of a. As an example, the A5 field completely decouples from the 4d 
physics also when one moves along the vertical line in Fig. [Jj (at a fixed /?s larger than the limiting 
value of the critical coupling which is reached for N5 — > 00), as was suggested in Refs. p~8 t fT9 | f56] . 
If one interprets the 5d model as an effective theory of a more fundamental theory defined with 
a cutoff a, this feature has a striking consequence: by ever increasing the extent of the extra- 
dimension in units of the cutoff, one ends up with a four-dimensional theory with no trace of its 
extra-dimensional ancestry. 

On the other hand, the cutoff can also be decoupled by approaching the phase transition line 
moving along a line of fixed m^/y/a (and therefore of fixed 7715, if one assumes that the 4d physics 
scale a is fixed): this is the case in which the resulting continuum theory features a scalar field 
coupled to 4d YM. To remove the cutoff in the 4d model at fixed m^/y/a one should study the 
lattice model at points along a line like, e.g., the brown curve in Fig. [JJ In this limit, going to 
larger and larger N§ values means that aa 2 tends to zero (continuum limit) and also oms goes to 
zero due to the approach of the second order transition line. However, the non-static Kaluza-Klein 
modes decouple in this limit, as it is only the static Kaluza-Klein mode which becomes critical and 
whose correlation length diverges at the second order transition. 

In short, our 5d lattice YM model with one compactificd dimension admits two different types 
of continuum limits in its dimensionally reduced phase: 

• ordinary 4d YM theory; 

• 4d YM theory coupled with an adjoint scalar field. 

Both limits are well-defined, renormalizable four-dimensional quantum field theories. 
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o p 5 

Figure 7: Sketch of the lines of constant low-energy physics (labelled by the m^/y/a ratio) in the 
(/?5, N$) phase diagram Fig. [5] of the 5d lattice model after taking the continuum limit in the fifth 
direction only. The red, leftmost, line indicates the phase transition separating the confining and 
the dimcnsionally reduced phases of the lattice theory. The lines of constant physics are displayed 
as solid black lines; at large j3§ values, they tend to hyperboles of equation P5N5 = const. The 
dashed blue lines correspond to fixed values of aa 2 : the lattice spacing a decreases when f3 5 —> 00 
at fixed N§ (green horizontal arrow), or when N§ — > 00 at fixed ,85 (green vertical arrow). In both 
these limits the A§ field decouples, and one recovers the ordinary continuum limit of the 4d YM 
theory On the contrary, a continuum limit including 4d YM coupled to an adjoint scalar Higgs 
can be obtained by moving upwards along a line of fixed m^/y/a (brown arrow). Finally, observe 
that, for any fixed, finite value of aa 2 (i.e. on any dashed blue line), there is a finite maximum 
value of A^5, i.e. a finite maximum value of L^/a, which is reached at the /?5 value where the dashed 
blue line hits the phase transition line. 
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6 Conclusions 



In this article, we have presented a non-perturbative lattice study of 5d SU(2) YM theory with 
one compact extra dimension. This model is a simple example in a class of theories which could 
be interesting for TcV physics and physics beyond the Standard Model, if large extra dimensions 
are realized in nature and can be probed by the forthcoming high-energy experiments at the Large 
Hadron Collider. 

At weak (bare) coupling and for all compactification lengths Ls/a, the model exhibits a pro- 
nounced scale separation, which grows exponentially as a function of L$/a, between the static 
KK- modes of the four-dimensional gauge fields and all the other modes. In this paper we have 
presented numerical evidence showing that this picture of dimensional reduction carries over to 
non-zero couplings as well. 

As the correlation length of the lightest modes is inevitably very large compared to the cutoff 
scale, the effects of the details of the regularization are strongly suppressed for low-energy observ- 
ables. In particular, the rotational symmetry is approximately restored. This defines a sensible 
scheme to obtain non-perturbative information about the low-energy sector from simulations of 
the five-dimensional lattice theory. 

In particular, using this approach it is possible to derive an upper bound for the radius of the 
extra-dimension for a given four-dimensional correlation length and ultra-violet cutoff (see Fig. [7j 
when era 2 is fixed, L^/a « N$ is bounded). For example, taking the electroweak scale as the 4d 
scale and assuming the cutoff to be at the Planck scale gives a maximum compactification length of 
the order of lOMp . Also, the large scale difference can potentially be used to acquire information 
about the compactification and cutoff scales from experimental measurements of four-dimensional 
quantities. If the mass 7715 and the self-coupling A of a scalar Higgs-like particle can be determined, 
say, at future LHC experiments, this information can be^ non-perturbatively mapped via lattice 
simulations to the two parameters of our lattice model ^5,^5, from which both the compactification 
scale and the cutoff scale can be estimated. For example, Fig. [5] illustrates how an experimental 
measurement of tos/ct 1 / 2 can be used as input to determine era 2 , and thereby the cutoff a -1 . 

Our main results are summarized in Fig.[2]and in Fig. [7] Fig.[2]shows the general structure of the 
phase diagram of the theory after extrapolation to the continuum limit only in the fifth direction, 
which is achieved by letting N$ — ¥ 00 (while keeping the 4d lattice spacing a finite), as a function 
of (3$ = 4a/.o| and 7V 5 = N 5 /j. Fig. [7] displays the behavior of the lines of constant low-energy 
physics (solid black curves) parametrized by constant values of the m§ / ^fa ratio (m^ being the 
non- vanishing mass associated with the static-modes of the field) in the dimensionally reduced 
phase. 

The cutoff can be sent to infinity in two different ways, both of which correspond to well- 
defined, renormalizable quantum field theories: 4d YM theory, and 4d YM theory coupled to an 
adjoint scalar field. The numerical results of our simulations are consistent with those of a similar 
study [2TJ[22], however our approach (based on a dynamical definition of the lattice spacing a 
from the four-dimensional string tension) leads to different physical conclusions. In particular, we 
emphasize that our findings are completely consistent with the non-renormalizability of the 5d YM 
theory One can construct continuum limits yielding ordinary four-dimensional models, while it 
is impossible to define a full-fledged continuum limit preserving the 5d nature of the YM theory, 
including the whole tower of Kaluza-Klein modefl 

In summary, the present study has clarified the phase structure of 5d YM theory, discussing the 
features of the phase diagram obtained from lattice simulations in a regime where a separation of 
scales makes the low-energy physics insensitive to the details of the theory at the cutoff scale. 

As for future lines of research, possible extensions of this study include: 

• extension to larger non-Abelian gauge groups; although it is reasonable to expect that the 
general qualitative features of 5d lattice YM theories are already captured by the present 

3 We have found no indication in our study for a non-trivial UV fixed point in a 5d isotropic system, as conjectured 
in |57H5t)l and investigated on the lattice in I16| 
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Figure 8: Finite-size scaling of the fourth-order Binder cumulant B4 = 1 — (P£}/(3(P§] 
N 4 x 2 lattice with 7 = 2. The horizontal line shows the mean- field value for the 4d Ising model 
and the vertical line indicates our determination of the critical value of (3. 



study of an SU(2) model, the extension to higher-rank groups might lead to a variety of 
spontaneous symmetry breaking patterns; 

• inclusion of fermionic degrees of freedom; in particular, a 5d model could provide a natural 
setup to study electroweak symmetry breaking and related problems. Furthermore, it would 
also be interesting to discuss the generalization of the setup studied in the present work by 
including lattice fermions in the domain- wall formulation |60[l61j. 
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Appendix: Finite-size scaling analysis of the phase transition 

To ascertain the order of the phase transition, we performed a finite-size scaling analysis for a 
iV 4 x 2 lattice at 7 = 2 for a range of volumes, measuring the fourth-order Binder cumulant 
associated with the average Polyakov loop along the fifth direction: 

The Binder cumulants measured on lattices with various volumes (shown in Fig. [8]) as a function 
of /3 are expected to cross at the critical point, at a universal value of B4 characteristic of the 
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Figure 9: On the left-hand panel, the reduced x 2 -value of the fitting ansatz Eq. (|A.5|) . applied to the 
data of Fig. [S] (the points corresponding to N — 4 were not included in the fit). On the right-hand 
panel, the data at j3 = 1.4826 are well described by the ansatz Eq. (|A.5[) , with x 2 /d.o.f = 0.50(7) 
and quickly deteriorates away from this point, which we take as our estimate of the critical point. 



universality class, provided that the transition is of second order. The transition breaks the Zi 
symmetry of the Polyakov loop and thus the universality class should be that of the 4d Ising model. 

Since four is the upper critical dimension for the Ising model, the critical value of the Binder 
cumulant can be estimated by considering a mean-field effective action with only one degree of 
freedom [52] : 

Z cS = [d<t> cxp(-S* cff ) (A.2) 

S e S 



N 4 Q^ 2 + ±u/) , (A.3) 



where t is the reduced temperature and 4> the mean field. 

Note that there is no contradiction between the effective action Eq. (|A.3[) where m / 0, and 
the triviality of </> 4 theory in four dimensions. The latter states that the renormalized coupling 
constant vanishes, whereas 4> in Eq. (|A.3[) stands for the bare field measured in the simulation: the 
mean P 5 Polyakov loop. 

At criticality, i.e. for t = 0, one gets: 

r 2 (i/4) 

Bt = 1 - - 0.27052. (A.4) 

Furthermore, the fact that four is the upper critical dimension for the Ising model also implies 
that the corrections to the critical value of B4 decrease only logarithmically with the linear extent 
of the system: 

Bi(N) = Bl + A/ln(7Y). (A.5) 

While rendering the direct observation of the crossings of the Binder cumulants numerically un- 
feasible, this form allows one to estimate a second-order critical point as the value of (3 where the 
iV-dcpcndence of the Binder cumulant is well described by the ansatz. Fig. [H] shows the reduced 
X 2 values from fits to Eq. (|A.5|) at various constant couplings /?, and the fit at the critical coupling. 
We conclude that our data are compatible with a second order phase transition in the universality 
class of the 4d Ising model at f3 c = 1.4826(1). 
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